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Abstract. We investigate the chiral extrapolation of the lattice data for the 
light-heavy meson hyperfine splittings D* — D and B* — B to the physical 
region for the light quark mass. The chiral loop corrections providing non- 
analytic behavior in m^r are consistent with chiral perturbation theory for 
heavy mesons. Since chiral loop corrections tend to decrease the already too 
low splittings obtained from linear extrapolation, we investigate two mod- 
els to guide the form of the analytic background behavior: the constituent 
quark potential model, and the covariant model of QCD based on the ladder- 
rainbow truncation of the Dyson-Schwinger equations. The extrapolated hy- 
perfine splittings remain clearly below the experimental values even allowing 
for the model dependence in the description of the analytic background. 



1 Introduction 

In the past few years there has been much progress in lattice gauge theory with 
many physical quantities having been calculated. Among them the hyperfine 
splittings in the heavy meson systems are of particular interest. With the aid of 
nonrelativistic QCD (NRQCD) on the lattice, the authors of Ref. ^ reported 
three lattice data for the qQ meson hyperfine splittings D* — D and B* — B. These 
data were obtained in the unphysical region where nig corresponds to being 
larger than about 680 MeV. With a naive linear extrapolation to the physical 
rriTj, the extrapolated hyperfine splittings are typically 120 MeV for D* — D, and 
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32 MeV for B* — B. These are significantly smaller than the experimental values 
140 MeV and 46 MeV respectively. The obvious shortcoming in the naive linear 
extrapolation is that the nonanalytic terms in the light quark mass, generated 
by chiral loops, do not appear. 

To correct this, we previously [21 IHj included pion loop contributions by ap- 
plying heavy meson chiral perturbation theory at small myr, following a series of 
works in this direction This leads to rapid, nonanalytic variation when 171-^ is 
smaller than about 400 - 500 MeV (which corresponds to a current quark mass 
niq ~ 60 MeV). When m^r is larger than 400 - 500 MeV, the heavy meson mass 
varies slowly and smoothly and is linear in m^, as indicated by the lattice data. 
Such considerations of the behavior of the individual meson masses led to the 
suggestion that the hyperfine splitting, yp = mp* — mp, for the heavy meson, 
P, can be extrapolated with the following form: 

yp = ap + ap + bpml, (1) 

where ap = ap* — ap, and ap(^pt^ the pion loop contribution to the meson mass 
of P (P*). This guarantees the correct chiral limit behavior, and ap and bp are 
fit parameters. The resulting extrapolated values for the hyperfine splittings for 
both D and B mesons are even smaller than those obtained in the naive linear 
extrapolation OIHj. 

Here we investigate the appropriateness of the last two terms of Eq. in 
representing the "analytic background" physics to which the chiral loop contribu- 
tions are added. We investigate the constraints from two models: the constituent 
quark potential model (CQM), and the Dyson-Schwinger equation (DSE) model 
in ladder-rainbow truncation "5". For all the hadron properties which have been 
calculated in lattice QCD, the lattice results vary slowly and smoothly when 
ruq is larger than about 60 MeV. This is characteristic of constituent quark be- 
havior [S] and suggests the use of the constituent quark mass for an efficient 
description of hadron properties in this region. The CQM emphasizes such a 
concept. 

On the other hand, the DSE ladder-rainbow model is a covariant modeling 
of QCD that satisfies the chiral symmetry constraints. In particular, it properly 
incorporates dynamical chiral symmetry breaking |7| which plays a dominant 
role in hyperfine splitting of physical ground states |Hl El- No explicit chiral 
loops are present at ladder-rainbow truncation; nevertheless the relation between 
pseudoscalar mass and quark current mass has the correct leading behavior at 
low and high quark mass 10 . 

An analytic background form for hyperfine splitting is deduced from the 
CQM in Sec. [5J and from the DSE model in Sec. and comparisons are made 
there. The chiral loop contributions to the hyperfine splitting are discussed in 
Sec. ^ and the chirally extrapolated results are presented in Sec. El A summary 
is made in Sec. El 

2 Constituent Quark Potential Model 

The constituent quark model (CQM) has been shown to work quite well even 
though it is a very simple model which has not been derived from Quantum 
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Chromodynamics. The constituent quark mass, Mq (where q denotes a, u ox d 
quark), is hnked to the current quark mass, m^, in the fohowing way j6j: 

Mq = M^ + c niq, (2) 

where c is of order 1 and is the constituent quark mass in the chiral hmit 
which originates from dynamical chiral symmetry breaking in QCD. 

In the constituent quark potential model, the potential between a quark and 
an antiquark consists of scalar and vector parts. The hyperfine splitting is caused 
by the relativistic spin-spin interaction between the quark and the antiquark. We 
consider the model where the vector part of the potential is caused by one-gluon- 
exchange. Then the hyperfine splitting for the heavy meson Qq is 

where Mq is the constituent quark mass of the heavy quark Q, is the strong 
coupling constant, and ip{0) is the wave function of the heavy meson at the origin. 
On the other hand, ip{0) can be expressed as IV'lO)]^ = ^ Mq+m (^^(^))- 
In a heavy meson which consists of a light quark and a heavy antiquark the dy- 
namics is mainly determined by the light degrees of freedom due to heavy quark 
symmetry. The long distance part of V{r) is the most relevant, and for a typical 
potential model with linear confinement, {V'{r)) is essentially independent of 
quark masses. Consequently, from the CQM we expect the mass dependence of 
hyperfine splitting to be well-represented by 

yCQM _ dp 

-Mq + MQ ' 

where dp is a constant. 

The constituent quark mass Mq in Eq. Q will vary with because of 
Eq. (jSJ and the variation of current mass niq with due to chiral symmetry. The 
latter variation we take from the Gell-Mann-Oakes-Renner (GMOR) relation in 

the form: = — ^fe-r, where rn^^^'' is the quark mass associated with the 

physical pion mass, m^^^^- Lattice studies indicate this relation is acceptable up 

phys 

to ~ 1 GeV . Then with Eq. a one has M„ = M^ + f^^^m^ 0, and 
Eq. @ takes the more convenient form 

y^^^ = ^^ . (5) 
ep + 

Here dp = dp ^ and ep = e + Mn phys with e being defined as 

c rUq c rriq 

('^Phys-|2 

n/f \^iT ) 

-'"Y phys • 

^ c rriq 

We expect Up^^ to be a better representation of the "analytic background" 
physics than the last two terms of Eq. for several reasons. Firstly, while 
^CQM g]^g^j,gg -^^^^j^ backgrouud part of Eq. Q a linear behavior with small 
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m^, its dependence on Mq is consistent with heavy quark effective theory [T^ . 
In contrast, Eq. ^ provides no exphcit Mq dependence. Secondly, the two pa- 
rameters of Eq. are totally phenomenological while the parameter ep of Up^^ 
is constrained by existing applications of the CQM to other baryon physics. For 
example, the ratio of the constituent quark masses of the strange quark and the 
u, d quark is given by = f_^^*'™hy'i° ^ , where {m\)s is the pion mass squared 

when ruq is equal to nis- Use of the CQM and chiral perturbation theory to guide 
the extrapolation of the magnetic moments of the spin- 1/2 baryon octet leads 
to Mg/Mu^d ~ 1.3 We let this ratio vary between 1.2 and 1.4. Use of the 
GMOR relation between and ruq leads to {m'^)s ~ 0.5 GeV^. Consequently 
we ahow the range of e' to be from 1.2 GeV^ to 2.5 GeV^. The CQM fit in 
Ref. j6, found c m^^^ ~ 5.90 MeV. If we take the heavy constituent masses to 
be Mc ~ 1.3 — 1.7 GeV and ~ 4.6 — 5.0 GeV, then we have the parameter 
ranges eo ^ 5.5 - 8.2 GeV^ and es ~ 16.5 - 19.2 GeV^ 

In Sec. 13 we add the pion loop contribution to yp^^ and fit the single free 
parameter dp to the lattice data, while constraining ep to be within the above 
expected range. For the sake of comparison, the dependence from alone 
is illustrated here by a preliminary fit of ?/p^^ to the lattice data, allowing both 
dp and ep to vary freely. This is displayed in Figs. Handle (The meaning of the 
model denoted CQM' in those figures is explained in the next Section.) We obtain 
do = 0.547 GeV^ en = 4.462 GeV^ and dp = 0.447 GeV^ ep = 14.71 GeV^ 
for y]^^ and y]^^ respectively. These ep values are just below the expected 
range. 



3 Dyson-Schvifinger Equation Model 

Information on hyperfine splitting in the absence of chiral loops can also be 
obtained from the modeling of mesons through the Dyson-Schwinger equations 
(DSE) of QCD in rainbow-ladder truncation |14L IH]. This a Poincare covariant 
approach that implements the correct one-loop renormalization group behavior of 
QCD 13 and the relevant Ward-Takahashi identities. The resulting axial current 
conservation preserves the Goldstone nature of the pseudoscalars jjj- With one 
infrared parameter (besides the quark current masses) the model ^1] provides 
an efficient description of the masses and decay constants of the light-quark 
pseudoscalar and vector mesons ^lE]) the elastic charge form factors Ft^{Q'^) 
and FxiQ"^) JSl and the electroweak transition form factors of the pseudoscalars 
and vectors jSl ^| • 

The current quark mass dependence of the DSE results satisfies the exact 
QCD pseudoscalar mass relation [Jj. For qQ mesons it takes the form j2I mp fp = 
{nig + uiq) Rp where the electroweak decay constant fp, and the projection Rp 
of the Bethe-Salpeter wave function onto 75 at the origin of qq separation, are 
dependent upon the current masses uig and tuq. The quark masses and Rp are 
dependent upon renormalization scale /i but the mass relation is not. The chiral 
limit of Rp is Rp^ = —{QCl)fi/fp, where fp is the chiral limit value; thus the 
GMOR relation follows as a collorary of the exact relation in the low mass limit. 
[The error in the GMOR relation at the K meson is 4%, at mq = 0.4 GeV the 
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Figure 1. The "analytic background" contribution to the D* — D mass splitting as obtained 
from the ladder-rainbow DSE model and from the CQM model. They are fitted to the lattice 
data for comparison purposes. 



error is 14%, while at the D meson the error is 30% 9 .] The heavy quark limit of 
the relation produces 10 mp oc uiq in accord with constituent quark behavior. 

The current quark mass dependence of the DSE results for the ground state 
nip and mp* has been studied recently 0121 with niq fixed at while uiq 
varies in the range 0-300 MeV. The results are summarized well in terms of a 
single parameter a by the forms 

nip = uq {nig + niq) + (m^ + niqf' (6) 
mp* = mp + a {niq + nig - 2mu) , (7) 

where current quark masses are defined at a scale of = 1 GeV. The large 
niq = mq behavior is consistent with recent DSE model studies jlOl I18j . Here 
ao = 1.724 GeV is the DSE model-exact value of —{qq)^/{fp)'^', use of this term 
alone corresponds to the GMOR relation. The term accounts for the cor- 
rection required by the exact pseudoscalar mass relation as reflected in the 
DSE results. The value a = 1.079 provides a good overall representation of the 
n/ p, K/K* , D/D* and B/B* meson masses. The absence of a linear term in 
Eq. respects dynamical chiral symmetry breaking and this in turn provides 
a good representation of the increase in the physical hyperfine splittings with 
decreasing mass. 

For any value of the parameter a, the hyperfine splitting obtained from 
Eqs. (jnj and in the form yp^^ = mp* — mp becomes independent of mg 
as mq oo. In that limit, heavy quark effective theory makes a more specific 
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statement, namely yp — > 0. This can be satisfied with the special value 



a 



4m,, 



{1 



1^ _ AmuOo 



(8) 



which yields a = 1.138. Since the employed DSE calculations are limited by 
iTiQ < 300 MeV, we prefer to use a = 1.138 as a way of incorporating large mg 
physics. The "analytic background" as deduced from the DSE model is thus fully 
proscribed in magnitude and m^r dependence. To make an initial assessment of 
this mTT-dependent background, in comparison to y^^^, we fit the form 



Up 



mp* — mp + Ap 



(9) 



to the lattice data with Ap being a free constant. The results are shown in Figs.^ 
andEl The required shifts Ad = -0.0244 GeV and Z\b = -0.0218 GeV perform 
a role which is later taken over by the mTr-dependent pion loop contributions. 
As we note later, the pion loop contribution in the vicinity of the lattice data is 
indeed negative and of this order. 
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Figure 2. The "analytic background" contribution to the B* — B mass splitting as obtained 
from the ladder-rainbow DSE model and from the CQM model. They are fitted to the lattice 
data for comparison purposes. 



We note that the conversion of mg-dependence into mjr-dependence, consis- 
tently within the DSE, has required the inversion of Eq. @ with mq = ruq. One 
may ask whether the resulting m^^^(m^) differs sufficiently from the GMOR 
result m^'^*-^^(m^) to warrant concern for the sensitive hyperfine splittings con- 
sidered here. We provide information on this in Figs. HandlSlby means of two 
versions of the y^^ description of the background. The version labelled CQM 
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employs m^^^^(m^); the version labelled CQM' employs m^^^(m^). The re- 
sults show this not to be a significant source of uncertaintity; it amounts to a 
maximum of only 3 MeV in yo at m.,^ = 0, and less for ys- Uncertainties in the 
DSE description of the background from the parameterization in Eqs. © and 
are not large. If we were to retain the parameter value a = 1.079 which is 
preferred by the low mass DSE calculations, the changes to the results shown in 
Figs, n and |2l are insignificant. 

From Figs. ^andl^Jit is evident that there is more curvature in yp^^ than 
in y^^ . This is due to the non-linearity imposed by chiral symmetry; the DSE 
model respects the axial vector Ward-Takahashi identity F71. In Tabled we sum- 
marize these representations of the "analytic background" in terms of would-be 
contributions to the hyperfine splittings at m-,^ = 0; we also include the linear 
extrapolation results. The DSE model and the CQM are consistent with each 
other to about 5% and both tend to be slightly below the naive linear extrapo- 
lation. However, the linearly extrapolated splittings are already significantly too 
low: 15% low in the D meson case and 30% low for the B meson. 



Table 1. Various background models extrapo- 
lated to m,r = and expressed as contributions to 
the hyperfine splitting. As explained in the text, 
m^^^^iml) is used in CQM, while m^^^{ml) is 
used in CQM'. All quantities are in GeV units. 





CQM 


CQM' 


DSE 


Linear 


Expt 


Vd 

tJB 


0.1224 
0.0304 


0.120 
0.0302 


0.1153 
0.0303 


0.120 
0.032 


0.140 
0.046 



4 Pion Loops 

The pion loop contributions to the masses of heavy vector and heavy pseudoscalar 
mesons were calculated in Ref. 0. The interaction between the pion and heavy 
mesons is generated from the following chiral Lagrangian for heavy mesons which 
is invariant under both chiral symmetry and heavy quark symmetry in the chiral 
and heavy quark limits, respectively jl9j : 

£ = -TT[Haiv^{D>'H)a] + gTT{HaH,j^Alj5), (10) 

where is the covariant derivative and is the axial-vector field, both of 
which include Goldstone boson fields, g is the coupling constant describing the 
interactions between heavy mesons and Goldstone bosons, and Ha (a=l, 2, 3 for 
u, d, s quarks, respectively) is a field operator which includes heavy pseudoscalar 

(P) and heavy vector (P*) mesons: Ha = -5^ (Pa ^7/^ — Pa 75) {v is the velocity of 
the heavy meson). The correction term ^-TiHacr^^ Ha(Tfj,p was added to Eq. (fHHl 
(where A2 is a parameter) since this term is responsible for hyperfine splittings. 
Based on experimental data we chose g"^ to vary between 0.3 and 0.5 and A2 
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Figure 3. Extrapolation of the lattice data for the D* — D hyperfine splitting using chiral 
loops, with a dipole form factor (solid line) and with a sharp-cutoff form factor (dashed line), 
added to the CQM model. The parameters are: A=OA GeV, A2 = -0.015 GeV^ = 0.4 and 
eD=6.9 GeV^ 



between -0.03 GeV^ and -0.02 GeV^ ^. With Eq. ^ ap can be expressed as 
the following: 

ap = - ^ „ / dfc- 



7r2/^ Jo ^/k^ + ml{y/k^ + ml - A) 



9' r^k.f^^'^Hk) 



4vr2/^ Jo + ml 

+ ^ (11) 
STT^f^Jo y^k^ + mli^k^ + ml + Ay 

where f-j^ = 132 MeV is the pion decay constant, A = —8X2/mQ, k is the 
absolute value of the three-momentum of the pion in the loop, and u{k) is an 
ultra-violet regulator. Since the leading nonanalytic contribution of pion loops 
is only associated with the infrared behavior of the integrals in Eq. (|11|) . it does 
not depend on the details of the regulator. In this work we choose two different 
forms for the regulator. One is the sharp-cutoff, 9{A — k), and the other the dipole 
form, which is more realistic, A^/^A^ + A;^)^. Here A characterizes the finite size 
of the source of the pion - i.e, the heavy meson's radius. In the fit we let A vary 
between 0.4 GeV and 0.6 GeV as in our previous work. 

The hyperfine splittings for D and B mesons were calculated in Ref. PP, 
where NRQCD was used to treat the heavy quarks. For the bare gauge coupling 
P = 5.7, the inverse lattice spacing is about 1.116 GeV. The box size is 2.1 
fm, corresponding to the volume 12'^ x 24. In their simulations, three different 
values for the hopping parameter k, 0.1380, 0.1390, and 0.1400, were used. The 
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light quark mass is related to k as niq = ^(1//^ — l/^c); with Kc = 0.1434. 
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0.0275 
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0.0265q q2 0.3 0.4 0.5 0.6 0.7 0.8 
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Figure 4. Extrapolation of the lattice data for the B* ~ B hyperfine splitting using chiral 
loops, with a dipole form factor (solid line) and with a sharp-cutoff form factor (dashed line), 
added to the CQM model. The parameters are the same as those in Fig.j^except that e_B = 17.9 
GeV^ 

The next step is to fit these lattice data with the form of Eq. (|12|) to determine 
the parameter d. The formulas in Eq. ()11() are obtained in the infinite volume 
limit. Since the lattice simulations are performed on a finite volume grid, the 
finite size effects should be taken into account. Following the second reference 
in Ref. |^ and Ref. [!20j we replace the continuum self-energy integral in Eq. 

by a sum over the discrete pion momenta which are allowed on the lattice: 
f d^k ss where the discrete momenta k^, ky, k^ are given by 

L is the number of lattice sites in each spatial direction, and the integer n satifies 
the constraint < n < ^. With 1/a = 1.116 GeV and L = 12, the smallest 
momentum allowed on the lattice, 2TT/aL, equals 0.58 GeV. 

With the least square fitting method we fix the parameter dp for different 
values of ep, A, g, and A2 in both the sharp-cutoff and the dipole schemes. It is 
noted that when < there is a pole in the first integral in Eq. (|11|) . In this 
case, we have kept the principal value of the integral which is real. The difference 
between the integral and its principal value is imaginary and corresponds to the 
width of the heavy meson. 

5 Results for Chiral Extrapolation 

Based on the above considerations, we extrapolate the lattice data with the form 




(12) 
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dp is a free parameter while, from Sec. |21 we constrain ep to be within e£) ~ 
5.5 - 8.2 GeV^ and eg ~ 16.5 - 19.2 GeV^. This direct extrapolation of the hy- 
perfine splitting should lead to smaller errors than the separate extrapolation of 
masses carried out in Ref. j^j . 

For ap, the allowed ranges of the parameters A, g, and A2 are those described 
in Sec ID In Figs. 01 and \^ we show the fits to the lattice data for A = 0.4 GeV 
and for the intermediate values of other parameters in both the dipole and the 
sharp-cutoff regulator schemes. When < 500 MeV the extrapolations start 
to deviate significantly from the constituent quark model. We require the fits to 
be within the uncertainties of the lattice data; this is always possible for A = 0.4 

GeV but not for A = 0.5 GeV or 0.6 GeV . In TableI21the hyperfine splittings yp, 
extrapolated to the physical pion mass, are listed as a range of values allowed 
by the variation of parameters as described above. We also give separately the 
relative uncertainties Ay / y caused by the errors in the lattice data. To assess the 
infiuence from the description of the background, we show results from use of the 
CQM background, Eq. (|12|). and from the linear background used in previous 
work, Eq. (^. The associated parameters are given in the Table. 

Table 2. The hyperfine splittings extrapolated to the physical pion mass through chiral loops 
added to the background described via the CQM (with parameter dp) or the linear repre- 
sentation (with parameters ap and bp). Results arising from the dipole form factor and the 
sharp-cutoff form factor are shown. The splittings yp are shown as a range of values allowed 
by the variation of parameters as discussed in the text. Separately given are the relative uncer- 
tainties Ay/y caused by the errors in the lattice data. 





D Mesons 


B Mesons 


Expt (GeV) 


0.140 


0.046 


Form Factor 


Dipole 


Sharp 


Dipole 


Sharp 


With CQM background 










yp (GeV) 
dp (GeV^) 
^y/y (%) 


0.0983-0.1150 
0.6675-1.0601 
1.24-1.45 


0.0978-0.1143 
0.6590-1.0237 
1.25-1.46 


0.0260-0.0289 
0.5048-0.6467 
1.82-2.02 


0.0259-0.0288 
0.4982-0.6245 
1.83-2.03 


With Linear background 










yp (GeV) 
ap (GeV) 
-bp (GeV-^) 
Ay/y (%) 


0.1000-0.1162 
0.1226-0.1395 
0.0308-0.0216 
6.28-7.30 


0.0965-0.1149 
0.1203-0.1375 
0.0340-0.0203 
6.36-7.57 


0.0270-0.0309 
0.0326-0.0374 
0.0077-0.0051 
8.90-10.18 


0.0256-0.0306 
0.0320-0.0369 
0.0086-0.0047 
9.00-10.76 



With our parameter ranges, the extrapolated hyperfine splitting with CQM 
background for D mesons varies from 0.0983 GeV to 0.115 GeV with the dipole 

'When the fit results are not required to be within the lattice uncertainties the high limiting 
values of yp remain unchanged while the low limiting values are reduced by about 8~17%. 
The values of the parameters dp, ap, and bp are not affected. 
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form factor, and from 0.0978 GeV to 0.114 GeV with the sharp-cutoff form factor. 
For B mesons, it varies from 0.0260 GeV to 0.0289 GeV and from 0.0259 GeV to 
0.0288 GeV, respectively. On the other hand, with use of the Hnear background 
the D meson result is in the range 0.100-0.116 GeV from the dipole form factor, 
and in the range 0.0965-0.115 GeV from the sharp-cutoff form factor. For B 
mesons, the extrapolated hyperfine splittings are 0.0270-0.0309 GeV and 0.0256- 
0.0306 GeV, respectively. The largest extrapolated hyperfine splittings we could 
obtain are about 115 MeV and 29 MeV for D and B mesons, respectively from 
the CQM background; with the linear background, these become 116 MeV and 31 
MeV, respectively. The present approach is unable to do better than 17% below 
experiment for the D splitting and 35% below experiment for the B splitting. 
These are further from experiment than the naive linear extrapolation which, 
as described in Sec. |21 produces a D splitting 15% below and a B splitting 30% 
below. 

One can compare the values at the midpoints of the ranges of the CQM 
background parameters ep and dp that fit the lattice data in the presence of the 
pion loop contribution, to the corresponding values obtained in Sec. [21 by the fit 
of the CQM background alone to the data. This confirms that at ~ 0.6 GeV^ 
the pion loop contribution to t/d is typically —10 MeV, and the contribution to 

is typically —20 MeV. These values are of the same order as the typical shifts 
Ap) and Ap found in Sec. |21to align the DSE background with the lattice data. 

In addition to the uncertainties which are associated with the parameters A2 , 
A, and ep we also analysed the uncertainties which are caused by the errors 
in the lattice data. As expected, these uncertainties are much smaller (less than 
2%) with the improved CQM background description compared to the linear 
treatment of the background (6~11%). Taking into account such uncertainties, 
we can see that the ranges of the extrapolated hyperfine splittings produced with 
the different background treatments are compatible with each other. 

Another interesting quantity is the slope of the fit line at large m-,^. From 
Figs. 01 and 21 we can see that the slope of the fit for B mesons is much smaller 
than that for D mesons. The CQM representation of the background carries 
a XjM.'q dependence for the slope and provides a natural fit. In contrast, the 
linear representation of the background contains no Mq dependence and this 
contributes to the larger uncertainties that result. The difference of slopes in 
these two extrapolation forms could be tested with more accurate lattice data. 

In the above fit we have treated the lattice data as truly statistical. Since 
the lattice data are highly correlated we could also fit the top and the bottom 
of the three data points in Figs. (31 and 01 to give the ranges of the fit. This may 
lead to a little change in the results. As an example, we consider the largest 
extrapolated hyperfine splittings for D and B mesons in Figures 01 and jH If we 
fit the top of the three lattice data in these two figures in the dipole scheme 
we obtain the extrapolated results as 0.1177 GeV and 0.0299 GeV for D and B 
mesons, respectively. These values are a little bigger than those obtained when 
we treat the lattice data as truly statistical. 
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6 Summary 

We have explored a number of issues that arise in the chiral extrapolation of the 
lattice data for hyperfine splittings of D and B mesons from the domain of large 
mass of the light quark to the physical point. We improve on previous work by in- 
vestigating the description of the non-chiral background via both the constituent 
quark model (CQM) and the modeling of QCD via the Dyson Schwinger equa- 
tions (DSEs) in ladder-rainbow truncation. We find that these descriptions of the 
background are consistent with each other. Thus the form employed here for the 
analytic background for extrapolation in is consistent with: chiral symmetry, 
heavy quark symmetry in the limit of infinite heavy quark mass, the constituent 
quark model, and the DSE modeling of QCD. We adopt the CQM background, 
add the difference of the chiral loop self-energies for vector and pseudoscalar 
mesons, and fit to the lattice data. We extract a range of extrapolated hyperfine 
splittings at the physical pion mass for both D and B mesons, as allowed by the 
theoretical uncertainties in parameters and pion-baryon form factors. 

In comparison to our previous work [21 |3j that employed an analytic back- 
ground linear in m^, the present approach has a better physical basis in the 
dependence of the background upon both Mq and m^. The pion-loop contribu- 
tion ap produces essentially all of the curvature evident in Figs. IHl and El at low 

and this is clearly less curvature than what we obtained in Ref. 2 . This is 
because the present approach deals directly with the hyperfine splitting and its 
extrapolation. There is a single form factor range A for a vector-pseudoscalar 
meson pair. In contrast, the prevous work |2l dealt with extrapolation of the me- 
son masses separately before the difference is formed. That procedure employed 
separate form factor parameters A for each meson and the resulting splittings 
were a difference of two quantities having different curvatures at low m^. The 
present more direct procedure has less curvature as a result of being better con- 
strained by physics. It has smaller uncertainties in the extrapolated splittings 
from errors in the lattice data. We conclude that the present procedure is more 
reliable. 

In contrast to extrapolations of lattice data for hadron masses j21j and mag- 
netic moments [221 ; the 1-loop chiral loop self-energies do not improve the agree- 
ment of the extrapolated results for hyperfine splittings in comparison with ex- 
periment. This work has explored the analytic background carefully and we con- 
clude that this is not a likely explanation. Although sub-leading nonanalytic 
behavior from higher order chiral loops could be examined, we conclude that 
more reliable lattice simulations for hyperfine splittings are required to resolve 
this situation. One would like to have data for full QCD, so as to not rely on 
the quenched approximation. The experimental mass splittings are small: 7% for 
D mesons and 0.8% for B mesons. The a ■ B term in NRQCD is of order 1/Mq 
and small corrections can have a magnified impact. The coefficient of this term 
may well be increased by the inclusion of radiative corrections beyond tadpole 
improvement, the possible light quark mass dependence of the clover coefficient 
in the clover action for light quarks, and higher order terms in NRQCD. With 
the physical improvements incorporated here, the extrapolation is unable to do 
better than 17% below experiment for the D splitting and 35% below experi- 
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ment for the B splitting. These are further from experiment than the naive hnear 
extrapolation which produces a D splitting 15% below and a B splitting 30% be- 
low. The fact that the deviation from experiment for B splitting is consistently 
double that for D splitting is suggestive of systematic errors in the lattice data. 
Higher quality lattice data would stimulate greater theoretical scrutiny of both 
analytic and nonanalytic behavior in for such a sensitive quantity. 
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